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1. INTRODUCTION 

The anyon [1] (quantum mechanical, non relativistic) model is a fascinating system to 
study. 

Firstly, it is a Aharonov-Bohm [2] system with no classical counterpart, in which planar 
particles interact at a distance topologically. Secondly, the particles being identical, the 
effect of the interaction is to make their statistics intermediate [1,3], nor Bose, neither 
Fermi. 

Little is known about anyons despite huge efforts devoted to their study. In the N- 
anyon case with harmonic attraction to the origin, linear eigenstates have been constructed 
[4], but they are known to be only part of the spectrum. Perturbative approaches [5] have 
undercovered a very complex structure for the equation of state. Finally, several numerical 
analysis have been developped [6]. 

Here, one is going to drop the indistinguishibility of the particles and propose a general 
model of particles topogically interacting at a distance a la Aharonov-Bohm. 

2. THE FORMALISM 

Consider the density of Lagrangian for N particles moving in plane minimally coupled 
to vector gauge fields A^{r) 

N ^ N N . 

Ln = Y. ^-^^^1 + E ^A^n)vi - Al{n))) + E ^ W y <9M^df (1) 

j=l OL=\ a,l3=l 

In the index i refers to the number affected to a given particle among N, the 
indices /x, i^, /), . . . correspond to the 3-dimensional Lorentz degrees freedom (/i = denotes 
the time direction, /i = 1, 2 the space directions), and a, /3, . . . label some internal degrees 
of freedom carried by the vector fields (1 < a < N). The Caj's are the electromagnetic 
couplings (topological charges) between the matter particles i and the gauge fields a at 
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position Tj. The Kap are the Chern-Simons^ self-couphngs of the gauge fields a,/3. 
One can easily see that under the gauge transformation 

= 9^A„ (2) 

the density of Lagrangian (||) changes by a time derivative, namely Je^jyp^J^S''^^ = 
d^i^fj.vp^ad'^ A^^)- Also, by the very definition (|l|) one can restrict to k^/j = Kpa, i.e 
to a symmetric matrix 

One proceeds by eliminating the time components of the gauge fields by varying 
the Lagrangian with respect to them [SL^ /5A°^ = 0). One gets 

TV N 

ep^5^{f- n) = J2 '^'apeoi^.d^A''^ (3) 

1=1 a=l 

where explicit use of the symmetry of [k] has been made. 

The magnetic field = eo^ud^A'^ appears in the right hand side of (^). Using the 
bidimensional identity d.^ = 27r5^(r) one finds in the Coulomb gauge d.Aa{r) = 

E-./3^>1 = E£^^X7^ (4) 

(k is the unit vector perpendicular to the plane). This equation can be symbollically rewrit- 
ten as [K][AQ,(r)] = [e/27r][A; x {f — fi) / (f — fi)'^] where [AQ,(r)] and [k x (r — rj)/(r — ri)^] 
are one column vectors on which the matrices [k] and [e] act. If one assumes that [k] is 
inversible one gets 

2tt [r — TiY 
The Hamiltonian corresponding to L^v stems from 



Pi = = mvi + V eaiAa{fi) (6) 

OVi „ 



^Here we concentrate on abelian Chern-Simons gauge fields. The generalization to the non-abehan case 
could be studied along the same lines. 
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One gets 

= ^Y.^Pi-H (^aiAc{fi)f (7) 
i a 

If one defines „4i(f) = Yl,a'^aiAa{r) one sees that = [e]*[K]~-^[e/27r][fcx(r — ri)/{r — rj)^] 

at position Vi enters the definition of Hm- 

Let us consider the coupUng matrix [a] = [e]*[K]~^ [e/27r]. It is by definition a symmetric 
matrix. The usual anyon model is nothing but taking [k\ and [e] to be onc-dimensional 

2 

matrices (i.e. one single gauge field) with a single anyonic coupling constant a = 
between the "fiux" = ^ and the "charge" e carried by each anyon. Singular self- 
interaction, which are present in A{f) at position r = fi, have to be left aside. 

Here, we get a general model where the couplings aij = aji can depend on i and j. 
One notes that Ai{r) is correctly defined at r = ri if and only if one asks for the matrix [a] 
to have its diagonal elements equal to 0. Thus one ends up with N(N — l)/2 independant 
anyonic coupling constant (the entries of [a]) and a Hamiltonian 

i=l 

~* A/ X V' ' 
where the gauge field Ai = Y^j^i c^ij with fij = fi — fj. 

In this more general point of view, the anyon model can be recovered by taking all the 
aij equal to a. But now, one gets as a bonus that .Aj(r)|^=^. is defined at the position of 
the particles, which was not the case in the original formulation. A quantum mechanical 
model of fiux tubes interacting a la Aharonov-Bohm with electric charges ej would 
correspond to taking e^j = implying (t)j = (Ea/jMa/s)^^ aij = ^ej(X;a/3Ma^)ej- 

[k] matrices have been already used [7] in order to reproduce fractionnal values for 
the Hall conductivity given as an = J2ai3['Aap- sghlightly different context (in 

particular there is no kinetic term for matter), one insists on a matrix [k] with integer 
entries in order to reproduce the quantum numbers (statistics, charge) of the electron. 
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Here the entries are not constrained to be integers, in [k] as well as in [e]. 

As in the anyon model, Ai^fi) is pure gauge, the singular gauge parameter being 



(9) 



k<l 



where O^i is the relative angle of the particles k and I. The strong analogy between the 
anyon model and the model proposed above will allow for the generalization of interesting 
results of the former to the latter. 

3. SOME EXACT AND PERTURBATIVE RESULTS 

i) Exact results 

The structure of the Hamiltonian HN{aij) given in ^ allows for the following general 
comment. Suppose one has an eigenstate ip{aij) of energy E{aij). By the virtue of the 
gauge transformation 



one finds that ip^aij) is a monovalued eigenstate of the Hamiltonian Hfyf{aij — niij) with 
the same energy E{aij). This implies that 'ip'{aij + mij) is a monovalued eigenstate of the 
original Hamiltonian H]\i{aij) with energy E{aij +mij). Thus as soon one knows an exact 
eigenstate, one can associate to it a tower of orbital eigenstates indexed by the quantum 
numbers rriij. 

Another consequence is that if one considers the m^j's as gauge parameter coefficients 
(not as actual quantum numbers), and choose rriij = E[aij\, where is the integer 

part of Uij, the gauge transformed Hamiltonian H]\j{aij — rriij) is then defined in terms 
of the couplings a[j = aij — E[aij\ with the constraint that a[j G [0,1]. Clearly, the 
physics described by both models is the same, thus one will always assume in the sequel 
that < aij < 1. In the anyon case with bosonic (fermionic) wavefunctions, the same 



ip'{aij) = exp{i^mkiOki)'tp{aij) 

k<l 



(10) 
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reasonning leads to < a < 2 (—1 < a < 1) since then the mjj's are constrained to be 
equal to a given even (odd) integer. 

The question is of finding particular eigenstates. Let us confine the particles by a 
harmonic attraction to the origin. This procedure [1,8] is commonly used in the anyon 
context since it yields a discrete spectrum. The A^-particle Hamiltonian with a harmonic 
interaction reads 

(11) 

One finds the relative eigenstates 

{M is a normalization factor, (3 = ^^j^ and = J2i<j ^ij) with eigenvalues 

{2n + N + -aij\)LO (13) 

i<j 

Note that since the states (12) are eigenstates of the total angular momentum operator, 
they are still eigenstates in the presence of a uniform magnetic field {uj — > y/uj'^ + lo"^). The 
integers ruij have to satisfy simultaneously either niij > (case I) or niij < (case 11)0. For 
a given N, the relative states < ri|n, mjj > have obviously too many quantum numbers. 
The TTT-ij-'s can be choosen as independant orbital quantum numbers, the other being either 
1 (case I) or (case II). One gets < n,mij\n' ,m'^j >= ^n,n'^m,m' where ra = J2i<j''^ijy 
and m' = J2i<j "^[j^ leading to sectors labelled by the quatum numbers n,m. In a given 
sector, the states can be separately orthonormalized. 

The states ( |l2|) with linear dependance on the ajj's narrow down to the usual linear 
anyonic eigenstates when one sets Uij = a. In this particular case one has additionnal 
conditions on the m^j's depending on the statistics (Bose or Fermi) imposed when a = 0. 
*If the Qij-'s are not constrained to be in the interval [0, 1], one gets either vriij > E[aij] or rriij < E[aij] 
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One way to reproduce (12) is to work [9] in the singular gauge 

ip' = exp(-i ^ akiOki)ip (14) 
k<l 

In complex coordinates = Xj + iyi the free gauge transformed Hamiltonian is 

/ 2 1 \ 

^'n = Y. [-—dA + -^muj'^z^Zij (15) 

The states 

i)' = exp(-^ ^ ZiZi)4> (16) 

i 

are eigenstates of the Hamiltonian (^) if (p is an homogeneous meromorphic function of 
degree d of zi, . . . , zn (case I) or of zi, . . . , zat (case II), with for eigenvalues {N+d)Lj. Since 
the prefactor exp(— J2i ^i^-i) can be rewritten as exp(— ^ ^i<j ^ij^ij) exp(— ^A^^ZZ), 
the center of mass coordinate Z = Zi/N factorizes out in (^) with energy lj. 

A suitable basis for (p is {JJi^j z^"^' °'^} or {Yli^^^j zfj^ where the m^j's are in- 

tegers, thus d = ^i^jifUij — Uij) or d = J2i<ji'^ij ~ ^^ij) (the m^j's are easily seen to 
be not independent since z™^'^ = {zij — zii )'"'•' can be expanded in a power series of zu). 
The space of eigenstates has to be a Hilbert space of square integrable functions where 
the Hamiltonian is self-adjoint. A simple requirement is to impose that the eigenstates 
have no divergence, implying that one has simultaneously rriij > (class I), or rriij < 



(class II). Thus one reproduces the eigenstates (12) with n = 0. An explicit calculation 
shows that non vanishing radial quantum number n > correspond to the Laguerre poly- 
nomials appearing in (|12|). The energy (2n -|- — 1 -|- J2i<ji''^ij ~ ciij))^ (class I) or 
i<j{aij — mij))uj (class II) coincides with ([T^). 
Imposing that the eigenstates vanish when any two particles come close together is 
an anyon-like requirment, amounting to the exclusion of the diagonal of the configuration 
space. Self-adjoint extensions [10] corresponding to diverging behavior at small distance 
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are possible. Indeed, allowing for short distance singularities implies that the m^j's have 
to satisfy simultaneously rriij > or rriij < 1. These states diverge at the origin but are 
still square integrable. Asking for the Hamiltonian to be self adjoint implies additionnal 
constraints leading to the possible self adjoint extansions rriij > or rriij ^ — li and 
rriij > 2 or rriij ^ 1- 

It is amusing to note that, in the situation where the a^j's are either equal to a or nul, 
the eigenstates (^) can be directly deduced from the anyonic eigenstates 

i<j 

These states are monovalued eigenstates of the N-anyon Hamiltonian if the integers rriij 
simultaneously satisfy rriij > E[a] or rriij ^ E[a] (here one cannot restrict a to be in the 
interval [0,2]). However if one drops the monovaluedness criterium one finds that these 
states are still solutions of the eigenvalue equation if some of the rriij are replaced by 
rriij = m'ij + Oi where the integers m'-- have to be simultaneously > or < 0. But one 
can get rid of the multivaluedness of the wavefunction by means of the singular gauge 
transformation exp(ia ^j^j where the last summation is performed only on those 
indices z, j for which rriij = Tn'ij + ot. One then reproduces the eigenstates defined above 
where the indices i, j for which Uij has been set to correspond to rriij = f^'ij + «• 
ii) Perturbative results 



Leaving aside the exact eigenstates (12), very little is known about the model defined 
above. In the case where the Ojj's are assumed to be small, a perturbative analysis can 
give some information on the system. Here again the experience gained in the study of 
the anyon model is helpful. A naive perturbative analysis might make no sense due to the 
very singulaij^ Aharonov-Bohm interaction afj/rfy One can circumvent this difficulty by 
^here we assume that the Hilbert space of unperturbed wavefunctions does not contain any states with 
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noticing that the singular gauge parameter Vt" = J2i<j '^ij^ij the imaginary part of the 
meromorphic function 

n{zi,Z2, ... ,zn) =^aij In Zij (18) 

i<j 

Let us "gauge transform" the Hamiltonian H]\f by taking as gauge parameter the real part 
of 0,, 0,' = J2i<j Oiij Inrjj i.e. 

V' = exp(±0')V^ = n rf^-'i'in, • • • , fTv) (19) 

i<j 

Because of the Cauchy-Riemann relations in 2 dimensions, ViQ" = k x Vjil' (implying 
dz^^ = iAk), one finds that the singular terms are absent in the Hamiltonian acting 
on ip 

Hn = E(|^ + E ^^^9. T E ^'4^^)■ (20) 
^ 2m ^ m rf, ^ m rf. 

As in the anyon model, one gets 2-body interactions with short distance behavior 
adapted to a perturbative analysis. One notes the it sign freedom in the choice of the 
redefinition of 'ijj. This sign freedom describes two possible short distance behaviors of the 
exact eigenstates as emphasized in the context of self adjoint extensions. The — sign in the 
redefinition (|T9|) corresponds to the self adjoint extansion rriij > (class I) or rriij < — 1 
(class II). 

iii) Some of the particles are fixed. 

So far the N particles are dynamical. It is however interesting to consider some of 
the particles fixed in the plane. To do so, take for the free Hamiltonian of the N-particle 
system in the singular gauge 

1 ^ 

^^'n = ^J:^^p' (21) 

i=l 



singular short distance behavior. 
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where Ci = 1 (moving particle) or = (fixed particle). In the Lagrange formulation 
(H), this amounts to affecting the speed vi with the factor e-j. Going back to the regular 



gauge via the gauge transformation (14) yields the desired Hamiltonian describing fixed 
and moving particles interacting via the Aharonov-Bohm couplings aij 

1 ^ 

HN = —y2e^{p,-Aif (22) 

1=1 

Most of the results presented above are still operative, in particular the perturbative 
considerations. However, the eigenstates ( p!2[ ) cannot be used anymore. 

Let us consider the scattering of a single particle by a finite lattice of — 1 identical 
flux tubes. This amounts to take ei = 1 and ej = for i = 2 . . . , N . Only the scattering 
by a single flux line is solvable [2]. In the case of two flux lines one has in the singular 
gauge 

tp' = exp{-iai20i2 - iai30i3)'ip (23) 

Let us assume that 2 and 3 are located at Z2 = —h and = h. For simplicity we 
set ai2 = ai3 = a- In the limit where the two flux lines are at the same point h ^ 0, 
one should reproduce the scattering of one electron by a flux line, where polar coordinates 
are used to separate the eigenstates equation. There exists a single coordinate system 
where the eigenstate's equation separates, and which contains the polar coordinates as a 
particular limiting case. This is the elliptic coordinates system [12]. Still working in the 
singular gauge but with the conformal mapping zi = /icosh(/x + where (^, (j)) are the 
elliptic coordinates, one gets the free Hamiltonian 

^' = ~2^h^cosh'' /X - cos2 0) + ^^^^ 



Notice that the gauge transformation (|23|) does not separate when one uses elliptic coor- 
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dinates. Indeed one gets 



/■/I -/IN sinh^ (fj. + i(f)) . „, 

exp(z^i2 + idis) = ^ (25) 

cosh /i — cos^ ^ 



This set of coordinates has been used in [11], however we stress that it does not 
describe the scattering of a charged particle by two isolated fixed flux lines, but instead 
the scattering of a charged particle by an elliptic flux tube. Indeed the eigenvalue equation 
is now separable and one can factorize the eigenstates as M{fj,)^((f)), leading to Mathieu's 
equations^ [12]. In the singular gauge, the angular function ^{(f>) is multivalued. With 
this choice of coordinates, the contour of the flux tube (which controls the multivaluedness 
of the singular free wavefunction) must necesseraly coincide with a geodesic defined by a 
constant /x. This is an ellips or in the singular case the line [—h, h] that connects Z2 = —h 
to ^3 = +h. It follows that the singular gauge transformation implied by the choice of 
elliptic coordinates is i/j' = exp(— iac/))^/'. 

Let us consider the singular flux [—h, h] line case. Mathieu's equations read 

- h'^mE cos{2(j))'^ + (A - h^mE)^ = (26) 

d^M + h^mE cosh(2^)M - (A - h^mE)M = (27) 

where A is a constant introduced to separate the coordinates. The general solution [12] 
for a multivalued $ is 

$ = Ae"0</' ^ C2r-e2*^<^ (28) 

r=— oo 

where iu^ is not an integer. A possible exp(— n^0) Yl,t^-oo C2r&~'^"'^ solution has been 
omitted since parity is broken anyway by the singular flux [— /i, /i] line. In the regular 
gauge, the wavefunction has to be monovalued when a complete winding encircling the 
®note that if a central harmonic interaction is added, Hill's equations have to be considered. 
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flux [—h,h] line is performed. This implies that (n,^ + ia)/i has to be an integer n. 
Obviously one can assume that n G [0, 2[. The quantification condition on A comes from 
the compatibility of the homogeneous system of linear algebraic equations determining the 
C2r's in terms of the coefficients 

mh'^E 1 



2 (2r - - A + mh'^ 

Thus one has two quantum numbers A and E, as desired (remember that has been 
fixed by the monovaluedness criterium). The solution of the equation on M introduces 
a parameter which does not yield any additionnal quantum number since there is 
accordingly a compatibility condition to satisfy^. 

In the limit where the flux [—h, h] line shrinks to a point, /i — > 0, one should reproduce 
the isolated flux tube case. One has /i — > oo, and zi — > re^'^ with r = {h/2) exp /i. The 
Mathieu's and compatibility equations become 

+ A$ = (29) 

- -drrdrM + = 2mEM (30) 
r r 

. 2 Oi-n . 2 ^ fo^^ 

sm vr — - — = sm (."jlj 

leading to the usual Bessel eigenfunctions with A = (a — £)^, where £ = n + 2p is the usual 
angular quantum number (remember that n = 0, 1). 

4. ON CURVED SPACE 

i) The formalism 

Let us consider a bidimensional manifold Ai2 defined by its metric gati^), where the 
^ one gets Ufj, = ±i(n — a). 
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indices a,b label the space coordinates The covariant Lagrangian is given by 

N 



1=1 \ a=l / o,/3=l 

(32) 

{g = I det {gab)\)- It is invariant under the gauge transformation 5A^^ = d^Aa- By definition 
[k] is a symmetric matrix. 

The Chern-Simons topological term having not explicit dependence on the metric one 
gets, by varying the Lagrangian with respect to A^^, an equation for the gauge potential 
identical to (^) 

N N 

ep,5\x - X,) = E '^"/^e^'^a^afe (33) 

1=1 a=l 

Defining Aia = X^a ^aiAaa the covariant Hamiltonian reads 

1 ^ 1 

^N = 1Z-Y. -^iPi<^ - ■^^a{Xi))g''''^/g{p^b - Aibix^)) (34) 

j^—i \/9 



where pia = jdia are the first quantized momenta (in the free case (34) is the Laplace- 
Beltrami operator). The gauge potentials Aia are determined by 

e''''daAb{x) =2nJ2 aij^\x - xj) (35) 
j 

with [a] = [e]*[K]~^ [e/27r]. The question is of solving this equation on a non trivial mani- 
fold. As on the plane one should proceed by first solving the 2-body problem 

e^'^diaAbixi, Xj) = 27rS^{xi - xj) (36) 

Then the vector potentials Aia{xi) are given as a sum of 2-body terms 

Aia{xi) = CXijAa{Xi, Xj^ (37) 
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plus possible irrotationnal one body terms describing the topology of the manifold (see 
below the cylinder case for an illustration). Aaixi,Xj) is the gradient of a potential 
symmetric under the exchange Xi — > xj. 

It follows that Aiaixi) is the gradient with respect to Xi of some function Q"{xi, . . . , xat) 
which is multivalued for any loop enclosing the singular points x = xj. In the singular 
gauge 

ip' = exp {-in") V' (38) 

the Hamiltonian is free but ip' is multivalued. 

In the Coulomb gauge da,/gA"' = 0, one can rewrite Af{xi) as the dual tensor 
— l/^e°'''dii)n'{xi, . . . ,xn). If one redefines 



tp = exp {±n') i> (39) 



the Hamiltonian acting on ip is 



^ / 1 1 1 1 \ 

= E o rP'ag^'VaPib + i-A'}{x^)dia T -^eabAnxi)d'^ (40) 

j^—i yzm^(^ mm j 

m and Q!' are by definition harmonic for the Laplace-Beltrami operator. In the simple 
case where Qah = ^^/gSab (plane, cylinder [13], sphere [14],...), Q' and satisfy the 
Cauchy-Riemann relation, and = + i^l" is meromorphic. 

ii) An example : the cylinder 

Let us consider the particular case of the cylinder R x Si . One has that and + 27r 
have to be identified. We define z = x^ + ix"^. Let us first construct Q. in the 2-body case 



(pq). On the one hand it must behave as aijinzij (see 18) when Zij — > 0, since locally 



the cylinder is equivalent to a plane. On the other hand dz^n = i{Aii — iAi2) has to be 
periodic in the variable xfj. One arrives at O = Oij ln(e^'J — 1) which, however, has yet to 
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be symmetrized under the exchange of i and j, yielding 



n = aij In 2 sinh ^ (41) 



One way to reproduce this result consists in considering the planar Jl' = aij In y {xjj)'^ + (xfj)'^ 
and making it periodic by introducing the infinite series 



oo 



n=— oo 



a,, Y: \nJ{x\f + {xl + 27:nY (42) 



This series is formally divergent, but it can be given a non ambiguous meaning by a usual 
procedure (derive with respect to x\j, perform the summation, and then integrate). One 
then obtains (p]). 

The essential difference between the plane and the cylinder consists in the non van- 
ishing contour integral of a gauge field on a non contractible loop around the cylinder. 
Consequently one has to introduce a line of flux inside the hole and, in the Coulomb 
gauge, consider a 1-body multivalued term (j)iZi. In the 2-body case this amounts to add 

(piZi + (pjZj to (H) 

Q, = aij In 2 sinh — - + (piZi + (j)jZj (43) 



Thus, in the A^-body case, one has 

_ ^ 

n= a*iln2sinh^^ — ^ + (44) 

0<i<j<N ^ i=l 

We now specialize to the 2-anyon case Q12 = a,(j)i = (j)2 = (j)- It follows that the 
relative motion is controlled by the anyonic a term, whereas the (j) term concerns the 
center of mass motion. 

Let us define the center of mass Z = {zi + Z2)/2 and relative z = zi — Z2 coordinates. 
In the singular gauge ( p8| ) the Hamiltonian is free. Since the configuration space of two 
identical particles is defined by the identification z ^ —z, the conformal mapping w = 
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2sinh(^/2) maps the cylinder on the plane and is thus well adapted. In polar coordinates 
w = re'^, the relative Hamiltonian reads 




Going back to the regular gauge amounts to the shift Oq ^ do — ia. 

Thus one has a situation similar to the relative scattering of two anyons on a plane. 
However, the Jacobian of the conformal mapping makes the relative motion non separable, 
implying seemingly out of reach exact cigenstates. 
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